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Abstract. The theory of first order Fermi acceleration at shocks 
assumes that particles diffuse due to scattering off slow-moving 
If} magnetic irregularities. However, cosmic rays are closely tied 
ON to magnetic field lines, and the transport process, particularly 
ON across the direction of the field, is likely to be more complicated. 
1—1 To describe cross field transport we employ recent extensions of 
Qnthe Rechester-Rosenbluth theory in localised stochastic regions 



of magnetic field. During acceleration at a shock and when 
the motion along the field is diffusive, there is a transition at 
1—1 a critical energy from "sub-diffusive" motion, where the mean 
T ~~ ' square displacement of a particle increases with time as t 1 / 2 , 
to compound diffusion, a combined process involving diffusion 
along a magnetic field which is itself wandering. Requiring this 
critical energy to be less than the system cut-off in a SNR of 
radius R places an upper limit on the coherence length, A, of 
Q the magnetic field for which diffusive shock acceleration can 
occur at quasi-perpendicular shocks when field line wandering 
ON dominates the effective transport of particles: A < R{U/v) x l 2 
4^ where U is the shock speed and v the particle speed. The accel- 
Q-ieration rate for particles in the sub-diffusive regime is derived 
q and its implication for the origin of high energy cosmic rays is 
S— i discussed. 

^ Key words: acceleration of particles - diffusion - plasmas - 
shock waves - (ISM:) cosmic rays - ISM: supernova remnants 



1. Introduction 

The theory of diffusive particle acceleration at shocks is based 
on multiple interactions of a particle with a shock front, between 
which the particle diffuses through the upstream or downstream 
plasma. In the simplest situations a power-law spectrum of ac- 
celerated particles is produced, independent of the diffusion 
coefficient. This process is thought to be responsible for the 
acceleration of cosmic rays to energies ~ 10 14 eV in super- 
nova remnants (see Blandford & Eichler 1987 for a review). 
The maximum energy which can be reached depends on the 
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particle's diffusion coefficient and the age of the shock. Most 
models of supernova remnants assume gyro-Bohm diffusion 
(usually referred to as simply "Bohm diffusion"), equivalent 
to large-angle elastic scattering with a mean-free path equal 
to the gyro-radius. Many shocks, such as those propagating in 
the wind of the progenitor star of a supernova explosion, are 
likely to be quasi-perpendicular, in which case the accelera- 
tion rate is determined by diffusion across the magnetic field 
(Ball & Kirk 1992). If the mean-free path is far greater than 
the gyro-radius then the spatial diffusion coefficient across the 
field lines k ± can be much smaller than the Bohm value, result- 
ing in a faster acceleration rate (Jokipii 1987). However, it is 
well-known that turbulence in a magnetised plasma can lead to 
wandering or braiding of the magnetic field lines: this has been 
extensively investigated in the context of transport in fusion 
plasmas but was first studied in connection with the propaga- 
tion of cosmic rays (Getmantsev 1963, Jokipii & Parker 1969). 
Clearly this process will influence the transport of charged par- 
ticles across the average field direction in the vicinity of a shock 
front (Achterberg & Ball 1994). Here we consider the situation 
in the light of recent work in transport theory (Isichenko 1991, 
Rax & White 1992, Chuvilgin & Ptuskin 1993) and find that 
there exists an energy range in which the transport of energetic 
particles cannot be described by a diffusion coefficient, but is 
rather "sub-diffusive" in that the mean square displacement of 
a particle perpendicular to the magnetic field increases with 
time as t x l 2 (cf. Getmantsev 1963). This has a strong influence 
on the acceleration rate, which we calculate in Sect. 4, using 
the propagator appropriate to sub-diffusion. At higher energies 
particle motion across the field becomes diffusive, the "com- 
pound diffusion" case, but with a coefficient which combines 
the effects of scattering off magnetic irregularities with mag- 
netic braiding. In Sect. 5 we discuss the implications of our 
results for the origin of high energy cosmic rays, and point out 
the main limitations of our treatment, which applies only in the 
quasi-linear regime. 

2. Braided magnetic field 

Consider a uniform magnetic field of strength B in the z direc- 
tion superimposed upon which there is a stationary stochastic 
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magnetic field characterised by a root mean square amplitude 
B\ and let b = B\/B. If the braiding of the field can be con- 
sidered quasi-linear, a field line passing through x = y = z = 
wanders such that its distance from the y-z plane x(s), as a 
function of the distance measured along the field line s (s = 
at z = 0), follows a Gaussian law: (a; 2 (s)) = 2Dms where 
(...) denotes an ensemble average, and Dm is a diffusion co- 
efficient, with the dimensions of a lengthscale. We consider 
turbulence characterised by the correlation lengths Ay and A_,_ 
along and perpendicular to the field, and by the Lyapunov length 
Al describing the exponential separation of adjacent field lines; 
i.e. field lines whose initial separation is less than A_,_. Then 
the diffusion coefficient of the magnetic field is D M = b 2 X l{ /4 
and the condition for quasi-linear behaviour is b <C Aj_/Ay 
(Kadomtsev & Pogutse 1979). Neglecting background shear of 
the magnetic field, quasi-linear theory gives A L = A^_/ (A N 6 2 ) 
(Isichenko 1991). 



3. Particle transport 

If the particles undergo large-angle elastic scatterings the re- 
sulting diffusion coefficients parallel and perpendicular to the 
field are related to the Bohm value k b = jv 2 mc/(3eB), for 
a particle of mass m, charge e moving at speed v (and with 
7 = (l-i» 2 /c 2 )" 1/2 )by 



K||=K B /e , K ± =€K B /(l + €) , 



(1) 



where e < 1 is the ratio of the collision rate to the gyro- 
frequency (e.g., Chuvilgin& Ptuskin 1993). In quasi-linear the- 
ory, the relationship between these quantities depends on the 
properties of the turbulence (e.g., Achatz et al. 1991). In par- 
ticular, e is in general energy dependent. However, it is usually 
assumed that both k ± and Ky are monotonically increasing func- 
tions of energy, and that k ± <C «h ■ In addition to this diffusion 
along and across the local magnetic field, a particle's position is 
also governed by the wandering or braiding of the field. Follow- 
ing Rechester & Rosenbluth (1978), consider a particle moving 
within a coherent patch of field lines with dimension A_,_ across 
the field. As the particle moves along the field, the patch is de- 
formed, but maintains constant area, provided the average field 
strength does not change. The field lines, however, diverge ex- 
ponentially on the scale of the Lyapunov length Al, so that the 
patch becomes filamented. After travelling a distance s along 
the field, the length of a patch filament is roughly A _,_ exp(s / Al), 
so that the typical thickness is d(s) = X ± exp(— s/Al). Thus, the 
patch is stretched and the filaments become so thin that the par- 
ticle can escape into an uncorrelated patch. Escape occurs either 
when the particle diffuses a distance d across the field, or when 
d becomes so small that the particle's gyromotion takes it out of 
the patch (Isichenko 1991). Thus, the escape or decorrelation 
time td 1S the solution of the equation 



s(td) = A L log 



Ax 



_Max(7- g , V^Kj.id). 



(2) 



We can now classify the type of transport according to whether 
the particle is able to travel this distance ballistically, or must 
diffuse along the field line. 

Ballistic motion, for which s = vt, requires the particle to 
propagate for a time given by the solution of Eq. (2) before its 
parallel velocity is affected by scatterings i.e., t& < t sc = k^/v 2 . 
According to Eq. (1) this condition implies r g > v^«I*d and 
hence vtd = Al log(X ± /r g ). In this case diffusive motion in 
the x direction starts after a time t = X n /v, and continues 
indefinitely. Thus, the particle suffers not only 'real' cross-field 
diffusion (for non-zero k ± ) but also transport due to field line 
wandering: 

(x 2 (t)) ( t X 2 \ 

= K± + D M v = K± I 1 + J^log^/rg) I (3) 

for = Al log(Xj_/r g )/v < t sc . This formula is valid for all 
t > X n /v and yields an effective diffusion coefficient which 
is large compared to k ± . The time at which the particle leaves 
a coherent patch of field plays no role, since diffusion is a 
Markov process i.e., it can be stopped and restarted in a new 
patch without any effect on the transport. 

If the particle remains within a coherent patch of field 
for a time t>K {l /v 2 , its motion in s is diffusive: s(t) = 
■^/2tK l{ . According to Eq. (1) this implies r g < ^4nj_td so that 
Eq. (2) yields y = log(A/y), where y = ^/2K||t d /A L and 
A = 6 2 A||/(eAj_\/2)- If A< 1, field-line wandering is unim- 
portant. There is then no anomalous transport regime, and 
the particle diffuses according to the coefficients of Eq. (1) 
(Isichenko 1991). Here, however, we are concerned with the 
opposite case. Assuming A > 1, the lowest order solution is 
y log A so that = A^log A) 2 /2/«|| which gives two pos- 
sibilities for motion perpendicular to the mean magnetic field. 
For times t < t d , the motion has a sub-diffusive component, 
since the particle is trapped within a coherent patch: 



(x 2 (t)) 



2t 



: K, + 



2 l/2 



D M K 



1/2 



tl/2 



for t < t<x 



(4) 



However, for t > the particle escapes and restarts motion 
in a new coherent patch. Over such timescales the transport 
process is Markovian, since memory of the previous patch is 
lost. The motion is diffusive: 



It 



2-PmK|| 
A l log A 



1 + 



log A 



(5) 



fort > t d = (A L logA) 2 /2/«|| > K\\/v 2 . This is the "collision- 
al" transport regime of Rechester & Rosenbluth (1978) and the 
process we refer to as "compound diffusion". 

4. Acceleration at shocks 

The existence of different transport regimes has been dis- 
cussed with reference to cosmic ray propagation, (Chuvilgin 
& Ptuskin 1993), where the behaviour on short timescales is 
not important. In the case of particles undergoing acceleration 
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at a shock front, however, these timescales must be compared 
with the average time spent by a particle in the upstream or 
downstream medium (the residence time, t les ) before return- 
ing to the shock front. For diffusive transport in the coordinate 
along the normal to the shock front t res sa K e ff/(Uv), where 
U is a fluid velocity and « e ff the effective diffusion coefficient 
(Toptyghin 1980). Thus, standard diffusion - governed by the 
coefficients of Eq. (1) - gives a residence time which is an 
increasing function of particle energy, leading to very slow ac- 
celeration at high energy. According to Eq. (5), the effective 
diffusion coefficient at a perpendicular shock in the compound 
diffusion regime is, apart from a logarithmic factor, proportional 
to K|| , which increases with energy. We expect this regime to ap- 
ply at sufficiently high particle energy, since the decorrelation 
time decreases with increasing energy whereas the residence 
time increases (see Fig. 1). 

Let us assume A^ = Ay = A, a macroscopic quantity com- 
parable with the scale of the accelerating system - e.g., the 
radius R of a supernova shock front. Eq. (1) assumes mag- 
netic irregularities of amplitude e sa ((6B) 2 } /B 2 on the scale 
of a gyro radius r g , which is much smaller than A (Chuvilgin 
& Ptuskin 1993). Transport is diffusive for those particles for 
which t res > t d , or, using Eq. (5) 

k ± A 2 A 2 (log A) 2 
Uv log A K|| 

for large A, where the diffusion coefficients n n and k ± are 
monotonically increasing functions of energy. This condition 
gives a critical energy below which the time a particle is trapped 
in a coherent field patch t& exceeds the residence time and the 
transport is sub-diffusive. The highest energy particles acceler- 
ated diffusively in a supernova can be found by setting the accel- 
eration time, t acc sa K e ff/U 2 , equal to the expansion timescale 
R/U. If acceleration to this energy occurs at predominantly 
perpendicular shocks where field line wandering is important, 
we arrive at the condition 

< * . (7) 
tf 2 logA U { ' 

From Eqs. (6) and (7) one can derive an upper limit 
fU\ 1 ' 2 R 

X < - n am/2 < 8 > 

\v J (log Ay/2 

which depends only logarithmically on b and e. If (8) is fulfilled, 
the highest energy particles can be accelerated diffusively. If it 
is violated, transport is sub-diffusive at all relevant energies. 
For energetic particles in supernova remnants U <C v, so this 
condition places a strong upper limit on the correlation length 
of the field fluctuations. 

A particle residing on either side of a perpendicular shock 
for a time t < t& will sub-diffuse when the cross-field transport 
is dominated by field line wandering. Such particles will never- 
theless be accelerated as they cross and re-cross the shock front. 
However, this problem has so far not been addressed, so that it is 



\ Sub — Diffusion 
\t d a p _1 


Diffusion 















Pcrit Pmax 



Fig. 1. Schematic view of the two regimes for particle acceleration at 
a quasiperpendicular shock. The decorrelation time is a decreasing 
function of particle momentum p so that particles with momentum 
below a certain value, p cr jt, sub-diffuse and diffusive shock acceler- 
ation is only appropriate for the highest energy particles for which 
Pcrit < p < Pmax- The case e = constant is shown. 

not a priori obvious that the standard results of diffusive shock 
acceleration (e.g., the spectrum and spatial dependence of the 
distribution) are unchanged. Here we restrict ourselves to a dis- 
cussion of the expected acceleration rate, which is different from 
that of the standard diffusive picture. In the steady state limit, 
with n s (p) the number density of particles at the shock with 
momentum p, the flux across the discontinuity is n s v/4 where 
v is the particle speed and the distribution is assumed isotropic. 
In diffusive shock acceleration there is a lengthscale, k,/U\, 
associated with the upstream distribution which is roughly the 
distance a particle can diffuse, x(t) sa \frd, before being over- 
taken by the shock, which is moving at speed U\ relative to 
the upstream medium. The total number of particles in the up- 
stream medium is then the number density at the shock times this 
lengthscale, ti s k/U\. Dividing by the flux gives t res = 4k/U\v 
as the mean time a particle resides upstream before it enters 
the downstream region. With a similar expression for particles 
in the downstream medium, moving at speed U2, which do not 
escape, the acceleration timescale for diffusive acceleration 

is obtained since a particle's momentum increases by a factor 
4(Z7] - U 2 )p/3v each cycle. In the sub-diffusive case, the dis- 
tance a particle can travel upstream in time t is now given by 
x(t) sa d\1 2 (/tyt) 1 ^ 4 . The lengthscale of the upstream distribu- 
tion, as defined above, is then D^[ 3 k^ 3 /U^ 3 . With the assump- 
tions of isotropy and time independence the mean residence 
time of a particle upstream is roughly (E^-Dm) 2 / 3 *;^ 3 jJJ\v to 
within a constant factor of order unity. Consequently, as far 
as the timescales are concerned and to order of magnitude, 
{U\Du) 2 I^Kyl 3 plays the role of an effective diffusivity for 
the non-Markovian transport, i.e. the acceleration timescale is 
roughly this expression divided by U 2 . While in the diffusive 

limit t res oc k„ , the residence time now increases less rapidly 

1/3 

with momentum, t KS oc n n ' (Fig. 1). 
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This argument can be made more formal by noting that 
the Gaussian propagator for diffusive processes is replaced by 
a probability density, P(x,t), of a particle being displaced a 
distance x in time t which depends on two processes. First there 
is the probability, Pf(x, s), that after a distance s along the field 
line the perpendicular displacement is x. Secondly P p (s, t) is the 
probability that scattering along the field transports the particle a 
distance s in time t. Then P(x, t) is the product of two Gaussians 
Pf(x, s) and P p (s, t) describing field line wandering and wave- 
induced diffusion, where s plays a temporal and spatial role 
respectively, integrated along the magnetic field line; 



P(x,t): 



d _ exp {-x 2 /4D M \s\) exp {-s 2 /4K n t) 



■M * 



The second moment of this expression gives (x 2 (t)) oc 

D\iK^ 2 t 1 / 2 as expected. Using the method of steepest descent 
the integral can be evaluated to give, after normalisation, 



P{x,t)K a(^D M Kl ,2 \x\t 1/2> j exp 



-/%| 4/3 



Z) 2/3 K 1/3 <l/3 



(ID 



(cf. Rax & White 1992) where a = 2" 1 / 3 (37r)- 1 / 2 and (3 = 
3/2 8 / 3 . The mean residence time of a particle upstream of a 
shock front is found by computing the total time t t spent up- 
stream of the plane x = U\t by a particle diffusing according to 
the propagator (10): t t = da; dt P(x, t). The flux escap- 
ing from upstream is vP{U\t, t)/4 which, when integrated over 
t > 0, gives the average number of times N CI a particle leaves 
the upstream region. The mean residence time is then 



^res — ^t/-^cr — 



4(U l D M ) 2 ' 3 Kl 
Uiv 



1/3 



(12) 



in accordance with the physical argument above. The accelera- 
tion timescale is then 



^acc — 



3D 2/ V /3 



1 



1 



u, 



1/3 



(13) 



with a trivial extension when D M and k„ are discontinuous 
across the shock. Eq. (13) only holds for particles which do not 
have time to decorrelate from the field in the vicinity of the 
shock, i.e. for the lower energy particles in Fig. 1 . If we assume 
A R, equating the expansion timescale of a supernova to the 
sub-diffusive acceleration time indicates that the highest energy 
to which particles can be accelerated is given by the solution 
of the implicit equation k b (RU)/(b 2 A) which replaces the 
usual equation valid for gyro-Bohm diffusion k b = RU. 

5. Discussion 

A long-standing problem in the theory of the origin of cosmic 
rays is that observations demand a single acceleration mech- 
anism for particles between lOGeV and 10 6 GeV. Diffusive 
shock acceleration in supernova remnants does not easily pro- 
duce 10 6 GeV particles in the time available. This problem is 



eased if supernovae accelerate particles at predominantly per- 
pendicular shocks, and if the particle transport at such a shock is 
governed by the diffusion coefficient k ± (Jokipii 1987). Such 
a picture requires an azimuthal magnetic field configuration 
around a supernova, which is perhaps not unrealistic (e.g., the 
spiral field swept out from a rotating star by a stellar wind). 
We have shown that if particles propagate ballistically over a 
decorrelation time the effective diffusion coefficient is large 
compared to k ± (Eq. 3). Furthermore, for particles which un- 
dergo scattering within a decorrelation time, we have shown that 
k ± is the relevant transport coefficient only if A < 1, i.e., only if 
the fluctuations on scales greater than r g are much smaller than 
those on this scale. Otherwise, field-line wandering controls the 
particle transport. The diffusive theory of particle acceleration 
can still be applied to the particles of highest energy, provided 
that the correlation length of the field fluctuations is sufficiently 
small. The transport of the highest energy particles is then com- 
pound diffusion, governed by an effective diffusion coefficient 
much larger than k ± . The problem of diffusively accelerating 
particles to 10 6 GeV thus remains. 

If the field is coherent on sufficiently long scales, the trans- 
port is sub-diffusive and a new theory of particle acceleration 
is needed. A preliminary calculation for the acceleration rate in 
this case indicates that rapid acceleration may be possible pro- 
vided b 2 A < 1 . However, the theory we have discussed has sev- 
eral limitations. We have neglected the effects of particle drifts 
and finite Alfven velocities on our estimates of the decorrela- 
tion time t&. These processes may aid a particle in decorrelating 
from a field fluctuation (Isichenko 1991) and thus make the 
regime of compound diffusion more accessible. We have also 
confined our considerations to the quasi-linear regime b <C 1, 
where it is possible to make simple phenomenological argu- 
ments on the lines of Rechester & Rosenbluth (1978) and have 
assumed an isotropic particle distribution at the shock front. It 
is unlikely that substantial progress can be made in relaxing 
these restrictions without undertaking numerical simulations. 
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